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(N ' Abstract. In this talk I present recent studies on vacuum polarization energies and 

& ' energy densities induced by QED flux tubes. I focus on comparing three and four 

dimensional scenarios and the discussion of various approximation schemes in view of 
the exact treatment. 
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(y-j . 1. Introduction and Motivation 
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In this talk I will present calculations of vacuum polarization energies that we pQ have 
performed for flux tube configurations in QED. Flux tubes in QED coupled to fermions 
exhibit interesting phenomena, such as the Aharonov-Bohm effect [2J, its consequences 
for fermion scattering |3] , parity anomalies [3] , formation of a condensate [5] , and exotic 
quantum numbers OIEHS]- Those (non-perturbative) features of the theory that give rise 
^| to these unusual phenomena make the analysis more difficult, especially in calculations 
that require renormalization. The investigation in ref. ^0] and the world line formalism 
in ref. [H] have addressed some of these issues. Here we provide a comprehensive 
approach employing techniques from scattering theory to analyze quantum energies of 
flux tubes. 

Our primary motivation for this analysis is to shed light on vortices in 
more complicated field theories, especially the Z-string in the standard electroweak 
theory |12j . The Z-string is a vortex configuration carrying magnetic flux in the field 
of the Z-gauge boson. Since the classical Z-string is known to be unstable ^3], the 
existence of such a vortex would require stabilization via quantum effects jHj, perhaps 
by trapping heavy quarks along the string. 

We compare the one-loop energies and energy densities of electromagnetic flux tubes 
in D = 2 + 1 and D = 3 + 1 spacetime dimensions. The classical calculation is of course 
the same in the two cases. The quantum corrections to the energy could possibly be 
very different ^T] because of the different divergence structure. In D = 3 + 1, the bare 
one-loop energy is divergent and requires renormalization. In D — 2 + 1, in contrast, 
the bare energy is finite. However, a comparison between the two dimensionalities is 
sensible only when we use the same renormalization conditions, which induces a finite 
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counterterm in the D = 2 + 1 case. Without this finite renormalization, the D = 2 + 1 
and D = 3 + 1 energies are qualitatively different. 

We also study this problem to analyze several technical puzzles associated with the 
computation of the one-loop energy of a flux tube. An efficient way to compute the 
energy is to use scattering data of fermions in the background of the flux tube. However, 
vortex configurations give long-range potentials, that do not satisfy standard conditions 
in scattering theory [T^], which usually guarantee the analytic properties of scattering 
data. In turn these properties are crucial to compute the vacuum polarization energy 
from scattering data. Hence we observe subtleties that emerge only because an isolated 
flux tube is unphysical, and once a region of return flux is included, the scattering 
problem is well-defined. In the limit where the return flux is infinitely spread out, the 
energy density becomes entirely localized at the original flux tube. 



2. Theory 

We consider the QED Lagrangian 

C = -\ (dpAv - d v A^f + $(ip- eA-m)iP , (1) 

where A^ is the Lorentz vector that represents the photon field and if) the spinor of a 
fermion which electric charge e. In D = 3 + 1 if) is required to have four components. 
Though that is not necessary in D = 2 + 1 we may nevertheless choose so. 
In radial gauge the flux tube configuration reads: 

* = 0, X B(r) = ^ d ^>, (2) 

where B denotes the magnetic field. For a GauBian flux tube we thus find 

fair) = 1 - e" r2/ " 2 , B G {r) = Baifte**^ , (3) 
which implies the net flux Fq = ttw 2 Bg{0) and the classical energy, 

Ed = \j d 2 rB 2 (f) (4) 

in D — 2 + 1. For D = 3 + 1 the identical expression gives the energy per unit length 
of the flux tube. We adopt the one loop approximation wherein the background fields 
consists only of gauge fields; so do the external lines of the Feynman diagrams. Thus 
we only consider fermions loops in the one loop approximation and the only relevant 
counterterm Lagrangian is 

C c t = ^V"- ( 5 ) 

It is natural to impose the renormalization condition that at zero momentum transfer the 
photon wave-function is not altered by quantum effects. This yields C^ 3 -* = — and 

tr 2 bnm 

Qifi^-e) — —j^-j (2 — 7 + In in D = 2 + 1 and D = 3 + 1 — e dimensions respectively. 
Note that the renormalization coefficient is finite in D = 2 + 1 while we have employed 
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dimensional regularization in D = 3 + 1. The corresponding counterterm energy (resp. 
energy per unit length) is 

= ^ [ d 2 xB> . (6) 



3. Fermion Determinant 

We obtain the vacuum polarization energy from the fermion loop by computing the 
functional determinants 

E<® = lim - [In Det(z^ - eA - m) - In Det(ip - m)] + E^ ] 

E<£.= lim -^-[lnDet(^-e4-m) -lnDet(^-m)] + . (7) 

Again, we consider the energy per unit length of the flux tube. Note that by inclusion of 
the counterterm contribution the above expressions are ultra-violet finite. In order to 
perform this computation we have to consider the non-trivial (static) background field 
Apix) 7^ of the flux tube, cf. eqs. (j3J) in the Dirac equation: 



a 



(p +A(x) S j + (5m m = uom. (8) 

The interaction induces a potential for the fluctuating fermions fields with two essential 
properties. First, bound states with energies Uj may emerge. In magnitude these 
energies are smaller than the fermion mass in*. Second, the continuum levels uj = 
±Vk 2 + m 2 acquire a non-zero phase shift, 8i(k) which translates into a modification 
of the level density, Ap(fc) = - J2e± ^ [<M^)]- Here i is the orbital angular momentum 
quantum number according to which the modes decouple and k is the linear momentum 
of the fluctuating field. We then find the vacuum polarization energy 

1 1 f°° k 

3 * 

Ef = -^Ef^ ln 4 + m2 -^ 2 )- A f^ln^EW 
d 8tt ^ \ J m? / 47r J m 2 v ' 

5~t(k) = 8i(k) — 5f\k) - 5f\k) . (9) 

In D = 3 + 1 the energy per unit length is obtained from the interface formalism 

of ref. It is important to stress that we have subtracted the first two orders 

of the Born series from the integrand to render the integrals finite. We will add 

back in these pieces in from of Feynman diagrams, E^. The identity between 

Born and the Feynman diagram contributions at a prescribed order has been verified 

within dimensional regularization ^H] and tested numerically, see e.g. appendix B of 

ref. jTSj. This procedure has the advantage that the renormalization conditions from 

For the flux tube configurations only threshold states with — m occur. The configurations that 
we consider in section 4 does not have any bound states. 
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the perturbative sector of the theory may be adopted [20j. To this end the renormalized 
vacuum polarization energy reads 

E[ D J = EP + E$8 + . (10) 



4. Embedding 

Configurations like eq. (jHJ) with non-zero net flux induce potentials that behave like 
V e ff(r) ~ 1/r 2 as r — > oo in the second order differential equations for the radial functions 
from which we extract the bound state energies and phase shifts. This behavior violates 
standard conditions necessary to deduce the analytic structure of scattering data in the 
complex momentum plane. As a direct consequence we observe that the phase shifts are 
discontinuous as k — > and Levinson's theorem cannot be employed to reliably predict 
the number of bound states. Though this is not a principle obstacle because we have 
other means to find the bound states and any singularity at k = is integrable in eq. (jOJ) 
it puts doubts on the use of scattering data for this computation. The analytic structure 
is furthermore mandatory to relate the matrix element of the energy momentum tensor 
to formulas like eqs. © and (fTUj) that underly our computation of the renormalized 
vacuum polarization energy. At the same time we observe that configurations with zero 
net flux are unrealistic. This becomes obvious from the Bianchi-identity: 

e Q ^<V^ = 0. (11) 

In D = 3 + 1 this identity comprises the well-known Maxwell equation d • B (x ) =0 
stating that magnetic field lines must be closed or extend to spatial infinity outside 
the region of interest. The latter scenario is not adequate for the study of the vacuum 
energy which requires to integrate over full space. In D = 2 + 1 the Bianchi identity 
becomes ^ = —d x E y +d y E x . This implies that it is impossible to create (static) net flux 
configurations from zero flux. This may cause inconsistencies as we want to compare 
energies of configurations with and without fluxes. We therefore superimpose a return 
flux configuration according to 

IQF G 



B R (r) 
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R 2 (1 + 256 (r 2 /R 2 - l) 2 ) (tt/2 + arctan(16)) 



B (r) = B G (r) + B R (r) . (12) 

In what follows we will refer to this configuration as the zero net flux configuration. 
It is straightforward to verify that lim^oo E c i[B Q ] = E C \[B G ] and lim^oo E-fd[B ] = 
Efb[Bg}- That is, the return flux does not contribute to the classical and counterterm 
energies as the position R of the return flux is sent to spatial infinity. The crucial 
question obviously is about the behavior of E vabC [B ] as R — » oo. To see what happens 
we compare the integrands of E$ for configurations with and without fluxes for two 
values of R in figure [T] The integrand corresponding to B oscillates around that 
corresponding to B G and these oscillations diminish as R increases. This indicates that 
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indeed lim^oo E vsx .[Bq] = E V&C [B G }. In order to unambiguously decide on that question 
we have to consider the energy density. This will allow us to distinguish between the 
contributions from the center flux tube and the return flux. 

5. Energy Density 

As motivated above, we consider the vacuum polarization energy density 

e(r) = 27rr(T 00 > (13) 

to decide whether or not a vacuum polarization energy can be attributed to a single flux 
tube. Here (T 00 ) denotes a specific matrix element of the energy momentum tensor T^ u 
in the background of the zero net flux configuration. As for the total energy the energy 
density is computed from three entries, 

e(r) = e*(r) + e FD (r) + e ct (r), (14) 

the contribution form scattering data, es(r) with appropriate Born terms subtracted to 
render the momentum integrals finite, the Feynman diagram contribution, epD( r ) as the 
Born subtractions must be added back in and the counterterm contribution from £ ct , 
eq. (jnj). For the details on this computation we refer to ref. pQ, in particular for the 
discussion on how the counterterm contribution cancels the UV divergences without the 
need for additional surface counterterms. Of course, a general consistency condition is 
that eq. (fTUj) is obtained from the spatial integral J °° dre{r). This requires to consider 
zero net flux configurations because only then the analytic properties of the scattering 
data are guaranteed that underly the proof of that identity, cf. ref. j2I] • 

In figure 121 we display the energy density as a function of the return flux position, R. 
A number of conclusions can be drawn from these numerical studies: 1) the energy 
density in central region of small r is independent of R, 2) the integrated density from 
that region gives lim^oo f^ 2 dr e(r) = E vac [Ba], 3) the energy from the return flux 
diminishes as its position is sent to infinity, i.e. lim^oo f^, 2 dr e(r) = 0. Altogether 
this is good news as it clearly confirms the naive method that only considers a single 
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Figure 1. The integrands of Ef \ eq. © with (black dashed lines) and without (red 
full lines) return fluxes for R = 6/m (left panel) and R — 26 /m (right panel). 
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central flux tube. The discontinuities in the phase shifts at k 
the vacuum polarization energy. 



do not propagate to 



6. Approximation Schemes for E vac [BG] 

Having established a reliable method for the computation of the vacuum polarization 
energies of flux tubes provides a good opportunity to employ this method to judge 
approximation schemes. Most popular are the derivative and perturbative expansion 
schemes. In both cases we will study the dependence of the vacuum polarization energy 
on the width w, introduced in eq. 0- It is important to consider variations of the 
background field that are consistent with the respective scheme. For example, for the 
derivative expansion to be appropriate we require configurations that vary slowly but 
keep the amplitude Bq(0) fixed. On the contrary, for the perturbative expansion we 
wish to consider different amplitudes and thus keep the magnetic flux Fq, fixed as we 
change w. The result for the derivative expansion are shown in figure 01 In both cases, 
D = 2 + 1 and D — 3 + 1, we find good agreement between the exact result and the 
leading contributions in the derivative expansion, even though the derivative expansion 
is known to be an asymptotic expansion and hence does not converge when summed to 
all orders jSj. In the perturbative expansion we evaluate the leading Feynman diagram, 

1 2 



E 



{D) 
FD 



8ttJ^ 
(4tt) d / 2 



dp 



°°, df(r) 
dr— — J (pr) 
dr 



dx 



x(l -x)pT(2-D/2) 
[m 2 + p 2 x(l — x)] 2 ~ D / 2 



(15) 



in D = 2 + 1 and D = 3 + 1 dimensions with the effective expansion parameter 
T = ^Fq- The UV divergence has not yet been removed from the diagram, eq. (|15j) 
and the counterterm part has to be added. This is illuminating especially for D = 2 + 1 
because with our renormalization condition it changes the leading behavior from T jw 2 
to T 2 /w A . Stated otherwise, renormalization is essential even for finite quantities. The 
change in this leading behavior is crucial to obtain agreement between the exact results 
and the perturbative expansion, both in D = 2 + 1 and D = 3 + 1, cf. figure |U 
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Figure 2. The left panel shows the energy density e(r) for various values of the return 
flux position R in D = 2 + 1. The right panel is the analog for the energy density per 
unit length of the vortex in D = 3 + 1. 
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Obviously the approximation schemes work well for the D = 2 + 1 and D = 3 + 1 
cases. Furthermore the vacuum polarization energies are negative and decrease with the 
width of the background flux tube for both cases. Also the energy densities (per unit 
length in D = 3+1) are similar in structure as can e.g. be seen from figure|21 As discussed 
above, to obtain these similarities it has been crucial to impose identical renormalization 
conditions. Having done so, it is to be expected that the simpler D = 2 + 1 case can be 
utilized to gain information about vacuum polarization energies of vortex configurations 
in more complicated D = 3 + 1 problems ^Bj • 




Figure 3. Renormalized one-loop energies in D = 2 + 1 (left panel) and D = 3 + 1 
(right panel), for fixed values of the magnetic field at the origin, as a function of the 
width of the Gaufiian flux tube. The lines correspond to the derivative expansion to 
second order The dots, circles, etc. represent the exact results for E vac . From top 
to bottom, eB G {0)/m 2 = 1.1,2,2.5,3,4,5. 




Figure 4. Renormalized fcrmion vacuum polarization energy in units of T 2 as a 
function of the width, for various fixed values of the flux T (2.5, 4.5, 6.5, 8.5 from top 
to bottom) in the Gaufiian flux tube. The full line represents the leading perturbation 
expansion contribution. The left panel is for D = 2+1 and the right panel for D = 3+1. 
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7. Conclusions 

In this talk I have reported on computations that we have performed for one-loop 
energies and energy densities of electromagnetic flux tubes in three and four spacetime 
dimensions. In general, this vacuum polarization energy contains ultraviolet divergences 
and an important feature of our approach is that it allows us to impose the standard 
renormalization conditions of perturbative quantum electrodynamics. Even though 
the calculation in three spacetime dimensions does not suffer from such divergences, 
a meaningful comparison between three and four dimensions can only be made when 
identical renormalization conditions are imposed. The use of scattering data to compute 
the vacuum polarization energy of an individual flux tube leads to subtleties arising 
from the long-range potential associated with the flux tube background, which does not 
satisfy the standard conditions of scattering theory. Consequently, the scattering data 
do not necessarily have the standard analytic properties required to relate the vacuum 
polarization energy to scattering data. We have therefore considered field configurations 
in which the flux tube is embedded with a well-separated return flux so that the total 
flux vanishes. We have constructed a limiting procedure in which this return flux does 
not contribute to the energy, enabling us to compute the energy of an isolated flux tube. 

We do not find qualitative differences between three and four dimensions for 
either the energy or energy density, once identical renormalization conditions have been 
imposed. However, we stress that renormalization in the case of three dimensions proved 
essential to this result because the (finite) counterterm contribution turned out to be 
large, thus causing sizable cancellations in the final result. 

This study gives an initial step toward understanding flux tubes and vortices in 
more complicated theories. In particular the similarities between the three and four 
dimensional cases can be used to determine whether quantum corrections stabilize the 
classically unstable strings in the standard model [16.. 
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